Super atom molecular orbitals: a new type of long-lived electronic states 
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We present ah initio calculations of the quasiparticle decay times in a Buckminsterfullerene based on the 
many-body perturbation theory. A particularly lucid representation arises when the broadening of the quasipar- 
ticle states is plotted in the angular momentum {t} and energy {e) coordinates. In this representation the main 
spectroscopic features of the fullerene consist of two occupied nearly parabolic bands, and delocalized plane- 
wave-like unoccupied states with a few long-lived electronic states (the superatom molecular orbitals, SAMOs) 
embedded in the continuum of Fermi-liquid states. SAMOs have been recently uncovered experimentally by 
M. Feng, J. Zhao, and H. Petek [Science 320, 359 (2008)] using scanning tunneling spectroscopy. The present 
calculations offer an explanation of their unusual stability and unveil their long-lived nature making them good 
candidates for applications in the molecular electronics. From the fundamental point of view these states illus- 
trate a concept of the Fock- space localization [B. L. Altshuler, Y. Gefen, A. Kamenev, and L. S. Levitov, Phys. 
Rev. Lett. 78, 2803 (1997)] with properties drastically different from the Fermi-liquid excitations. 
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Superatom molecular orbitals (SAlMOs) were recently dis- 
covered^ as universal characteristics of Ceo molecules and 
their aggregates. Further experiments on endohedral systems^ 
and calculations for a series of quasi- spherical molecules^ 
showed that SAlMOs are associated with the whole system 
(rather then with a particular atom), have a well defined spher- 
ical symmetry (Fig.[T] angular momentum i = 0, 1, 2, prin- 
cipal quantum number n = 3) and are capable of forming 
chemical bonds. Being markedly different from other unoc- 
cupied delocalized states they hold a promise for unique ap- 
plications in molecular electronics. Hitherto theory was not 
able to answer why these states are so resilient to the chem- 
ical environment and why their character is not washed out 
by the hybridization: questions of critical importance for the 
SAlMO-mediated charge transport^. 

To resolve these issues we report here results of quantum 
chemical calculations for prototypical Ceo using recently de- 
veloped self-energy formalism^ . When a transport electron is 
injected into the unoccupied SAIMO (the single-particle (Ip) 
state a) the molecule undergoes a transition into the mixed 
quantum state that decays in time. At later instances the de- 
cay is exponential exp(— 7^t) with 7^ given by the imagi- 
nary part of the on-shell self-energy 7^ = ^I]Q,cK(e) (Ref.[6]). 
However, it is the initial stage of the decay that is manifested 
in ultrafast nonequilibrium processes of technological rele- 
vance. We have shown^ that the spectral function for a large 
class of relevant electronic systems can be represented in the 
form 

A„(t) = exp(-7„^), (1) 

where the set-in time of the exponential decay is given by 
Ta = "^la/^aa- ^hc spcctral function Aa{t) has the fol- 
lowing short and long time-limits: 

f^AUt) n° -alj, (2) 
A„(t)*^-e-7.t, (3) 
meaning that for short times (t <C Tq,) we have the quadratic 
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FIG. 1. Single-particle states of Ceo molecule. Each group of the 
degenerate states is represented by a cylinder and is characterized 
by the averaged orbital momentum (x-axis), the energy (^/-axis), the 
energy uncertainty a (atomic units (a.u.), z-axis or colour coding) 
and the orbital momentum uncertainty (Si, the area of the cylinder 
base). The occupied valence states (as calculated fully ab-initio) 
form two nearly parabolic bands that can be described analytically 
by En{i) = En{0) + n = 1, 2. This is a signature 

of the delocalized motion of the electrons on the spherical cage of 
the radius R = 3.57 A (Ref.^. SAMOs show comparable values 
of and much smaller a and form a third band dissolving in the 
sea of other unoccupied states. Note that the weak dispersion of the 
SAMOs with i, as observed above, is due to their larger averaged 
radius as compared to the first two angular subbands. 
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FIG. 2. Diagrams for the electron self-energy: The exact expression 
in terms of the 2plh response function (G^^^^) in which the entrance 
and the exit channels cannot be separated by cutting one fermion line 
can be factorized approximately into the exact particle-hole (p-h) 
propagator (H) and the exact Ip Green's function (G) or, equiva- 
lently, into a product of the exact screened Coulomb interaction (W) 
and the exact Green's function. 

decay 

Aa{t)^l-alj^/2. (4) 

o-^^ is the central quantity for our theory. It can be de- 
termined as follows. A simple substitution of asymptotic ex- 
pansions into the Dyson equation leads to the exact relations 
between the spectral moments of the self-energy T,{uj) and of 
the spectral function A(a;) (Ref.|7l): 

M^^) =1, = M^^) - s, (5) 

5^(0)^M^2)_j^(i)j2^ (6) 

where Sloo is the frequency independent real part of the self- 
energy^, s is a diagonal matrix with the elements given by the 



where pij = (0|c|cj|0) is the single-particle density matrix 
and the Coulomb matrix elements are defined as (a/3 1 7(5) = 
/ (i(rir2)^a(ri)^/3(ri)(/)^(r2)(/)5(r2)/|ri -rsl- The prefac- 
tor 2 arises from the spin degeneracy. Apart from the fac- 
torization of the six-point correlation function Eq. ^ is ex- 
act. It is sufficient to determine only the density matrix of the 
correlated ground state. For small molecules or clusters this 
can be done at any desired level of accuracy by the config- 
uration interaction (CI) approach which, however, currently 
not feasible for fullerenes. Therefore, we consider a further 
approximation widely known as GW^ (Ref . TT) and treat the 
screened Coulomb interaction on the random phase approxi- 
mation level. This is equivalent to the configuration interac- 
tion singles (CIS) treatment of the excited states. In view of 
the Brillouin's theorem^"^ the ground state density matrix is 
Pij = riiSij leading to the energy-uncertainty expressible in 
terms of the Coulomb matrix elements only: 

i j,k 

The main computational burden is the transformation of the 
Coulomb matrix elements from the atomic to molecular or- 
bital basis which scales as 0{N^) in the brute force imple- 



zeroth-order state energies. By defining the matrix of spectral 
functions in terms of the imaginary part of the single-particle 
Green function (A{u;) = ^|ImG(co')|) and likewise for the 
spectral function of the self-energy (S(co') = |Im5](co')|), 
and by the use of the superconvergence theorem^ one can re- 
define the matrices in terms of the frequency integrals: 

/oo 
dujuj''A{uj), n = 0...2, (7) 
-oo 

/CO 
duiS{ui). (8) 
-CXD 

Equation ^ expressed in the basis of Hartree-Fock states 
is central for our discussion. In order to be feasible it should 
be written in a finite basis and coupled with an appropriate 
approximation for the self-energy operator. It is well known 
that the self-energy can be recast in terms of the 2plh (2hlp) 
six-point correlation function^ Our approach corresponds 
to its factorization into the product of p-h and p(h) correlation 
fun ction s (Fig. [2]) linking it to the established GW approxima- 
tioiPni. It is remarkable that after all complicated Lehmann- 
type expressions for the correlation functions are substituted 
into the T, = iGW and the frequency integral in Eq. ^ is 
performed we obtain an expression for the energy-uncertainty 
of a particularly simple form: 
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mentation. To obtain convergent results we used very large 
basis sets (up to 6-311++G(3df,3pd)) resulting in TV = 2340 
functions for Ceo- Integral transformations for this number 
of basis functions is not feasible with any standard quantum 
chemistry package and required a parallelized implementation 
fully accounting for the symmetry of the system. 

Our calculations clearly indicate a peculiarity of SAMOs: 
a strong localization in the energy domain (cf. tTHOMO = 
16.0 eV and trsAMOa =6.1 eV) (Fig. [T]) or their extended 
life-time even in comparison with the life-time of HOMO or 
LUMO states. This finding endorses the potential of SAMOs 
as transport channels in molecular electronic devices, since 
the energy is hardly dissipated during the short transport time. 
Our results are also unexpected from the Landau's theory of 
Fermi liquids and illustrate that finite systems possess elec- 
tronic excitations that differ drastically from quasiparticles in 
extended matter. Namely, a particle excitation may decay into 
two particles and one hole. This process may recur for many 
generations or it may stop after a few. For the former case 
a large number of many-particle peaks form a Lorentzian en- 
velope, the so-called quasiparticle. If the decay stops after a 
finite number of branchings, only Dirac peaks appear in the 
single particle spectrum (denoting a localization in the Fock 
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= Mg]-[M(i)]^^ = 5]2] = 2^ ^ {ai\jk){pi\j'k')[njnj>5kj6k'j' ^nk{l-n^^^^^^ (9) 

i j,k,k,k' 
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space). A decisive discriminating factor was shown to be the 
particle's energy (e) (Ref. 16 ). Our ab initio approach suggests 
that it is actually the kinetic energy ex that governs the decay. 
We find that even though the energies of electrons with differ- 
ent radial character may be very similar their kinetic energy 
is substantially different (as follows from the specific form of 
the Kohn-Sham potential for this system, Ref. 17 ), which im- 
poses strong restrictions on the Coulomb matrix elements and 



leads to the localization of SAMOs in the Fock space. This 
hints on the prolonged life-times and stipulates the observed 
high stability of these states. 
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